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We show using covariant techniques that the Einstein static universe containing a perfect fluid is 
always neutrally stable against small inhomogeneous vector and tensor perturbations and neutrally 
stable against adiabatic scalar density inhomogeneities so long as Cg > i, and unstable otherwise. 
We also show that the stability is not significantly changed by the presence of a self-interacting 
scalar field source, but we find that spatially homogeneous Bianchi type IX modes destabilise an 
Einstein static universe. The implications of these results for the initial state of the universe and its 
pre-inflationary evolution are also discussed. 



I. INTRODUCTION 

The possibility that the universe might have started 
out in an asymptotically Einstein static state has recently 
been considered in [1], thus reviving the Eddington- 
Lemaitre cosmology, but this time within the inflationary 
universe context. It is therefore useful to investigate sta- 
bility of the family of Einstein static universes. In 1930, 
Eddington [2] showed instability against spatially homo- 
geneous and isotropic perturbations, and since then the 
Einstein static model has been widely considered to be 
unstable to gravitational collapse or expansion. Never- 
theless, the later work of Harrison and Gibbons on the 
entropy and the stability of this universe reveals that the 
issue is not as clear-cut as Newtonian intuition suggests. 

In 1967, Harrison [3] showed that all physical inhomo- 
geneous modes are oscillatory in a radiation-filled Ein- 
stein static model. Later, Gibbons [4] showed stability of 
a fluid-filled Einstein static model against conformal met- 
ric perturbations, provided that the sound speed satisfies 
Cs = dp/ dp > i. 

The compactness of the Einstein static universe, with 
an associated maximum wavelength, is at the root of 
this "non-Newtonian" stability: the Jean's length is a 
significant fraction of the maximum scale, and the maxi- 
mum scale itself is greater than the largest physical wave- 
length. Here we generalize Gibbons' results to include 
general scalar, vector, and tensor perturbations and a 
self-interacting scalar field. We also show that the Ein- 
stein static model is unstable to spatially homogeneous 
gravitational wave perturbations within the Bianchi type 
IX class of spatially homogeneous universes. 

Consider a Friedmann universe containing a scalar field 
(j) with energy density and pressure given by 

and a perfect fiuid with energy density p and pressure p = 
wp, where — i < w < 1. The cosmological constant A is 



absorbed into the potential V. The fluid has a barotropic 
equation of state p — p{p), with sound speed given by 
Cg — dp I dp. The total equation of state is Wt — Pt/ Pt = 

{P + P<p)/{P + P<p)- 

Assuming no interactions between fiuid and field, 
they separately obey the energy conservation and Klein- 
Gordon equations. 



p + 3{l + w)Hp = 0, 
The Raychaudhuri field equation 
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(1 + 3w)p + </>2 - V{(t)) 



has the Friedmann equation as a first integral. 
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where K = Q, ±1, and together they imply 
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II. THE EINSTEIN STATIC UNIVERSE 

The Einstein static universe is characterized by if = 1 , 
d = a = . It is usually viewed as a fiuid model with 
a cosmological constant that is given a priori as a fixed 
universal constant, and this is the view taken in previous 
stability investigations [2-4]. However, in the context of 
infiationary cosmology, where the Einstein static model 
may be seen as an initial state, we require a scalar field 
(j) as well as a fiuid, and A is not an a priori constant 
but determined by the vacuum energy of the scalar field. 
The vacuum energy is determined in turn by the potential 
V{(j)) of the scalar field, which is given by some physical 
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model. An initial Einstein static state of the universe 
arises if the field starts out in an equilibrium position: 



r(0o)=O, A = 87rGF(0o)- 



(6) 



We will first discuss the simple case where V is trivial, 
i.e., a flat potential. A general Einstein static model has 
p = = V = (j) hy Eqs. (l)-(5), and satisfies 



-{l + 3w)po + 4'l = Vo, 
(1 + w)po + 4>l 
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(7) 
(8) 



If the field kinetic energy vanishes, that is if (/iq = 0, then 
by Eq. (8), (1 + w)pQ > 0, so that there must be fluid in 
order to keep the universe static; but if the static universe 
has only a scalar field, that is if po = 0, then the field 
must have nonzero (but constant) kinetic energy [5,1]. 
Qualitatively, this means that there must be effective to- 
tal kinetic energy, i.e. Pt + Pt > 0, in order to balance 
the curvature energy in the absence of expansion or 
contraction, as shown by Eq. (5). 
Equations (7) and (8) imply 



Wt = 
W4, = 



Vo + {l + 3w)po/2' 



Vo-il + 3w)po/6, 



(9) 
(10) 

and Po ^ 



It follows that < since w > 
0. In other words, the presence of matter with (non- 
inflationary) pressure drives the equation of state of the 
scalar field below — |. If the field has no kinetic energy, 
then by Eqs. (7) and (8), its equation of state is that of 
a cosmological constant. 



60 = w^ = -1. 



(11) 



If there is no fluid (po = 0), then the field has maximal 
kinetic energy. 



Po = 



(12) 



In this case, the field rolls at constant speed along a flat 

potential, V{(j)) = Vq. Dynamically, this pure scalar field 
case is equivalent to the case oiw — I pure fluid with cos- 
mological constant. This can be seen as follows: equating 
the radii of the pure-field and fluid-plus-A cases, using 
Eq. (8), implies Vo = {1 + w)po, by Eq. (7). This then 
leads to = -(3 + 5w)/(5 + 3w), by Eq. (10). How- 
ever, we know that = —1 for the fluid-plus-A case, 
and equating the two forms of leads to w = 1. 

Thus there is a simple one-component type of Einstein 
static model, which may be realised by perfect fluid uni- 
verses with cosmological constant (the pure scalar field 
case is cqiiivalcnt to the w; = 1 fluid model). This is 
the case considered by Eddington, Harrison and Gib- 
bons. The general case, of most relevance to cosmology. 



involves a scalar field with nontrivial potential V{(j)), so 
that the Einstein static model is an initial state, cor- 
responding to an equilibrium position {Vq = 0). This 
general case is a two-component model, since in addition 
to the scalar field, a fiuid is necessary to provide kinetic 
energy and keep the initial model static. 



III. INHOMOGENEOUS PERTURBATIONS OF 
THE FLUID MODEL 

Wc consider first the effect of inhomogencous density 
perturbations on the simple one-component fluid mod- 
els. Density perturbations of a general Priedmann uni- 
verse arc described in the l-|-3-covariant gauge-invariant 
approach by A = a^D^p/p, where D^ is the covariant 
spatial Laplacian. The evolution of A is given by [6] 



A+{2-6w + 3cl) HA + 
+ A-kG (3u;^ -I- Qcl - 8w - 
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12 [w - cl) 
l)p+ {2,cl - 5w)A] A 



- c^^D^A - «; ( D^ 3^ ) f = , 
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(13) 



where £ is the entropy perturbation, defined for a one- 
component source by 



p£ ^a^D'^p- pc^A. 



(14) 



For a simple (one-component perfect-fluid) Einstein 
static background model, f = and Eq. (13) reduces 
to 



Ak = 47rG(l + w) [1 + (3 - eyl] poAfe , 



(15) 



where we have decomposed into Fourier modes with co- 
moving index k (so that ^ — fc^/ao). It follows that 
Aft is oscillating and non-growing, i.e. the Einstein static 
universe is stable against gravitational collapse, if and 
only if 



(fc^-3)c^>l 



(16) 



In spatially closed universes the spectrum of modes is dis- 
crete, fc^ = n(n + 2), where the comoving wavenumbcr n 
takes the values n = 1, 2, 3 ... in order that the harmon- 
ics be single- valued on the sphere [7,3]. Formally, n = 
gives a spatially homogeneous mode (A: = 0), correspond- 
ing to a change in the background, and the violation of 
the stability condition in Eq. (16) is consistent with the 
known result [2] that the Einstein static models collapse 
or expand under such perturbations. This follows from 
the Raychaudhuri equation (3). 

For the first inhomogencous mode (n = 1), the sta- 
bility condition is also violated, so that the universe is 
unstable on the corresponding scale if there is initial 
data satisfying the constraint equations. However this 
is a gauge mode which is ruled out by the constraint 
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equations {Goi = SirGToi). This mode reflects a free- 
dom to change the 4-velocity of fundamental observers. 
(For multiple fluids with different velocities, one can in 
principle have isocurvature modes with n = 1 that are 
not forced to zero by the Einstein constraint equations.) 
The physical modes thus have n > 2. Stability of all 
these modes is guaranteed if 



(17) 



This condition was also found by Gibbons [4] in the re- 
stricted case of conformal metric perturbations; we have 
generalized the result to arbitrary adiabatic density per- 
turbations. Harrison [3] demonstrated stability for a 
radiation-filled model and instability of the dust-filled 
model, which are included in the above result. 

Thus an Einstein static universe with a fluid that satis- 
fies Eq. (17) (and with no scalar field), is neutrally stable 
against adiabatic density perturbations of the fluid for all 
allowed inhomogeneous modes. The case w = 1 covers the 
pure scalar field (no accompanying fiuid) models. When 
Cg > the stability of the model is guaranteed by Eq. (16) 
for all but a finite number of modes. The universe be- 
comes increasingly unstable as the fluid pressure drops. 
Clearly, a dust Einstein static universe (c^ = 0) is always 
unstable. In that case, Eq. (15) reads = AnGpoAk, 
implying exponential growth of (the Jeans instability; 
see [3]). 

The physical explanation of this rather unexpected 
stability lies in the Jeans length associated with the 
model [4]. Although there are always unstable modes 
(i.e., with wavelength above the Jeans scale) in a flat 
space, in a closed universe there is an upper limit on 
wavelength. It turns out that, for sufficiently large speed 
of sound, all physical wavelengths fall below the Jeans 
length. By Eq. (8), the maximum wavelength 27rao de- 
pends on the equation of state and is given by 



Gpoil + w) 



(18) 



The Jeans length is 27rao/nj, where A^. = and fc? = 
nj(nj -|- 2). Hence, by Eq. (15), we have 



(19) 



Stability means A < Aj, which leads to Eq. (17). For dust, 
Aj = and all the modes are clearly unstable. For radi- 
ation, Aj = 0.61Amax, and for a stiff fluid Aj = 0.81Am£ix- 
In both of these cases the Jeans length comprises a con- 
siderable portion of the size of the universe, and is greater 
than all allowed wavelengths. This would also place im- 
portant restrictions on the evolution of non- linear density 
inhomogeneities by shock damping or black hole forma- 
tion. 

We also note that the stability condition entails neu- 
tral stability; the oscillations in A are not damped by 



expansion, since the background is static. In the cos- 
mological context, where an initial Einstein static state 
begins to expand under a homogeneous perturbation, the 
expansion will damp the inhomogeneous perturbations. 
However, the Einstein static will not be an attractor; in- 
stead, the attractor will be de Sitter spacetime. 

Vector perturbations of a fluid are governed by the co- 
moving dimensionless vorticity Wa = acOa, whose modes 
satisfy the propagation equation 



= - {l- 3Cg) Hzuk . 



(20) 



For a fluid Einstein static background, this reduces to 



(21) 



so that any initial vector perturbations remain frozen. 
Thus there is neutral stability against vector perturbations 
for all equations of state on all scales. 

Gravitational-wave perturbations of a perfect fluid 
may be described in the covariant approach [8] by the 
comoving dimensionless transverse-traceless shear 'Eab = 
(i(^abi whose modes satisfy 



Sfe + SHUk + 
SwG 



- + 2- 



^ (l+3«;)p+-A 



Sfc = 0. 



(22) 



For the Einstein static background, this becomes 

Efe+47rGpo(fc' + 2)(l + w;)Efe =0, (23) 

so that there is neutral stability against tensor perturba- 
tions for all equations of state on all scales. However this 
analysis does not cover spatially homogeneous modes. 
It turns out that there are various unstable spatially 
homogeneous anisotropic modes, for example a Bianchi 
type IX mode which we discuss below. The associated 
anisotropies can die away in the case where the insta- 
bility results in expansion. If this is the case, despite 
the extra spatially homogeneous unstable modes, the ex- 
panding inflationary universe will be an attractor. 

In summary, for the simplest models (one-component 
perturbations), we find neutral stability on all physical 
inhomogeneous scales against adiabatic density pertur- 
bations if Cg > i, and against vector and tensor per- 
turbations for any Cg and w, thus generalizing previous 
results. 



IV. SCALAR-FIELD PERTURBATIONS 

We turn now to the case of a self-interacting scalar 
field, i.e., where the scalar field is governed by a non- 
fiat potential V{(j)), with initial Einstein static state at 
<f) — (po, i-e., Vq = = 4'o- This is the general dynami- 
cal problem: the stability analysis of the Einstein static 
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universe as an initial equilibrium position within a phys- 
ically motivated (non-flat) potential. (Some realizations 
of this scenario are discussed in [1].) Since we are inter- 
ested only in the behaviour close to the Einstein static 
solution, we can treat H and V as small in the pertur- 
bation equations. The lowest-order solution for density 

perturbations, A^°\ thus corresponds to i? = = V, 
and is given by the fluid perturbation solution, Eq. (15), 



lS.^^\t) ^ Akco&ujot + sin Wot , 
= 47rG(l + w)pQ [{e - 3)c^ - 1] 



(24) 
(25) 



with Ak and Bk constants for each mode. The next or- 
der contains contributions from the scalar field perturba- 
tions, 



(1) 



.(1) 



(26) 



during the time when the background is close to the Ein- 
stein static equilibrium position. 

The entropy term in Eq. (13) has no intrinsic fluid 
contribution since we assume the fluid is adiabatic. A 
scalar field generically has intrinsic entropy perturba- 
tions, which follow from Eq. (14) as [9] 



A.. 



(27) 



(See [10] for the corresponding expression in the metric- 
based perturbation formalism.) These entropy perturba- 
tions have a stabilizing effect on density perturbations of 
the scalar field: the entropy term in Eq. (13) cancels the 
preceding Laplacian term — c^D^A^, which would pro- 
duce instability when c| < 0; what remains is the term 
— D^A0, which contributes to stability. 

For an initial Einstein static state, the entropy contri- 
bution from the scalar field to lowest order is. 



(1) 

k 



2A 



(1) 



(28) 



where we used [c^]^"^ = — 1. There is also a relative en- 
tropy contribution, arising from the relative velocity be- 
tween the field and the fluid. We expect this contribution 
to be negligible. 

We now expand all the dynamical quantities in 
Eq. (13) to incorporate the lowest order effect of the self- 
interacting scalar field. In particular, a = ag + a(i) and 

p = Po+ where = ■^(p^^-^ + V(i), and similarly for 
the pressure. Then t« — > w -|- , where 



(29) 



and Cg ^ Cg -|- c^^^, where c^^^ is determined by the form 
of the potential V((}>) near (pQ. It also follows from the 
background equations that 



SttG 



2-^(1) +%) 



(1) = -V('i) . 



(30) 

(31) 



Using these results, Eq. (13) gives the evolution equation 
for the contribution from scalar field perturbations: 



A 



(1) 



(1) 
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(0) 



where 



= 4^G(1 + w)po [{k^ - 3)(c2 + 2) - 1] 
Fit) = -i2-6w + 3cl)H^i), 

G{t) = 1 {l2[c^,) - + 24[w - c^]^ 



(32) 

(33) 
(34) 



-|-6ww(i) + 6c(i) - 8«;(i) 
+ [3u;^ - Su; -I- 6c^ - 1] 



Po 



'(1) 



2c: 



0(1) 

ao 



(35) 



The change in the frequency, loq — s- ujq, shows the stabi- 
lizing effect of scalar-field entropy perturbations, which 
effectively increase the adiabatic sound speed term in ujq: 
Cg ^ Cg + 2. The solution of Eq. (32) is given by Green's 

method in the form 



(1) cos (Dot 



Wo 

sin LJot 

Wo 



j A^°^ [{FsinQoty - Gsinwo^] dt 
j a[."^ [{FcosCjQty -Gcoswot] dt. (36) 



Since Cjq > 0, for fc^ > 3, we see that stability is not 
changed by the introduction of scalar-field perturbations. 
We note that if general relativity is extended to include 
higher-order curvature corrections to the gravitational 
Lagrangian, then the existence and stability conditions 
for the Einstein static miiversc arc changed in interest- 
ing ways [12] which are linked to the situation of general 
relativity plus a pure scalar field through the conformal 
equivalence of the two problems [13]. 



V. SPATIALLY HOMOGENEOUS TENSOR 
PERTURBATIONS 

We now investigate the stability of the Einstein 
static universe to spatially homogeneous gravitational- 
wave perturbations of the Bianchi type IX kind. The 
anisotropy in these spatially homogeneous perturbations, 
which is absent within the Priedmann family, is what al- 
lows for tensor modes. 

The Einstein static universe is a particular exact so- 
lution of the Bianchi type IX, or Mixmaster, universe 
containing a perfect fluid and a cosmological constant. 
The Mixmaster is a spatially homogeneous closed (com- 
pact space sections) universe of the most general type. It 
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contains the closed isotropic Fricdmann universes as par- 
ticular subcases when a fluid is present. Physically, the 
Mixmaster universe arises from the addition of expansion 
anisotropy and 3-curvature anisotropy to the Friedmann 
universe. It displays chaotic behaviour on approach to 
the initial and final singularities if w < 1. This is closely 
linked to the fact that, despite being a closed universe, 
its spatial 3-curvature is negative except when it is close 
to isotropy. 

The diagonal type IX universe has three expansion 
scale factors ai{t), determined by the Einstein equations, 
which are [11]: 



((210203)' ^ 4[(a| - aly - gf] 
010203 (010.203)^ 

+ A + 47rG(l - w)p, 



(37) 



and the two equations obtained by the cyclic interchanges 
Oi — > 02 — > 03 ^ Oi , together with the constraint 

^ + ^ + ^=A_47rG(l+3«;)p, (38) 

Ol 02 03 

and the perfect fluid conservation equation (with w con- 
stant) 



p(0i0203) 



w+l 



const . 



(39) 



The Einstein static model is the particudar solution with 
all o, = oo = const. We consider the stability of this 
static isotropic solution by linearizing the type IX equa- 
tions about it: ai{t) = ao + Sai{t) , p{t) = po + 5p{t), and 
6p = w6p. 

The linearized field equations lead to 



12 8 
Stti H — 5- ^Oj = -5- (Sai + Sa2 + Sa^) 
do Oo 

-|-47rG(l - 'w)Sp. 

The conservation equation to linear order gives 

A \ /(5ai -|- (5o2 -|- ^03^ 



(5/9= -4 



SttG 



■ Po 



Oo 



(40) 



(41) 



If we define an arithmetic-mean perturbed scale factor 
by 



SA{t) = 



Sui + (5o2 + ^03 

Oo 



(42) 



then it obeys 



5 A = 35 A 



-(l-w)(A + 87rGpo) 



^0 



"0 



■3m;). 



(43) 



The first term on the right-hand side arises from the 
pure spatially homogeneous gravitational wave modes of 
Bianchi type IX (i.e., with 6p = = dp). The second 



term on the right-hand side arises from the matter per- 
turbations. Thus we see that the perturbation grows and 
the Einstein static solution is unstable to spatially homo- 
geneous pure gravitational- wave perturbations of Bianchi 
type IX, 



5 A oc exp ( 2t 



(44) 



When the matter perturbations are included the insta- 
bility remains unless 1 -I- 3w < 0, which we have ruled 
out for a fluid. In this case, the perturbations oscillate. 
This condition corresponds to a violation of the strong 
energy condition and this ensures that a Mixmaster uni- 
verse containing perfect fluid matter will expand forever 
and approach isotropy. Notice that the matter effects 
disappear at this order when w = 1. This is a famil- 
iar situation in anisotropic cosmologies where a w = 1 
fluid behaves on average like a simple form of anisotropy 
"energy" . 

The instability to spatially homogeneous gravitational 
wave modes is not surprising. Mixmaster perturba- 
tions allow small distortions of the Einstein static solu- 
tion to occur which conserve the volume but distort the 
shape. Some directions expand whilst others contract. 
Note that these S'0(3)-invariant homogeneous anisotropy 
modes are different to the inhomogeneous modes consid- 
ered in the perturbation analysis. Mixmaster oscillatory 
behaviour is not picked out by the eigenfunction expan- 
sions of perturbations of Friedmann models. 

The relationship between the stability of the Ein- 
stein static to inhomogeneous and spatially homogeneous 
gravitational-wave perturbations can be considered in 
the light of the relationships between the Bianchi type 
modes and inhomogeneous modes. In open universes, 
Lukash [14] has pointed out the correspondence between 
Bianchi type VII anisotropy modes and the inhomo- 
geneous gravitational wave perturbation spectrum that 
emerges when appropriate eignfunctions are chosen for 
solutions of the Helmholz equation on negatively curved 
spaces. These Bianchi type modes correspond to choos- 
ing complex wavenumbcrs [15,16]. In the case of a closed 
universe of Bianchi type IX a similar characterisation of 
the homogeneous tensor mode as arising by choice of 
an imaginary wave number would lead to the Mixmas- 
ter instability found above in the case of k'^ < —2. If 
such modes are admitted in the spectrum of perturba- 
tion modes for the closed geometry, then they will lead 
to instability. However, since superpositions of them lead 
to non self-adjointncss of the Laplacian [17], there will be 
problems with quantum analogues of these modes and it 
is not clear that they are physically admissible in the real 
universe. 
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VI. CONCLUSIONS 

There is considerable interest in the existence of pre- 
ferred initial states for the universe and in the existence 
of stationary cosmological models. So far this interest 
has focussed almost entirely upon the de Sitter universe 
as a possible initial state, future attractor, or global sta- 
tionary state for an eternal inflationary universe. Of the 
other two homogeneous spacetimes, the Einstein static 
provides an interesting candidate to explore whether it 
could play any role in the past evolution of our Universe. 
It is important to know whether it can provide a natural 
initial state for a past eternal universe, whether it allows 
the universe to evolve away from this state, and whether 
under any circumstances it can act as an attractor for 
the very early evolution of the universe. We might also 
ask whether it is not possible for it to provide the glob- 
ally static background state for an inhomogeneous eternal 
universe in which local regions undergo expansion or con- 
traction, manifesting an instability of the Einstein static 
universe. With these questions in mind we have investi- 
gated in detail the situations under which Einstein static 
universe is stable and unstable. 

We have shown that the Einstein static universe is neu- 
trally stable against inhomogeneous vector and tensor 
linear perturbations, and against scalar density pertur- 
bations if Cg > |, extending earlier results of Gibbons 
for purely conformal density perturbations. However, we 
find that spatially homogeneous gravitational-wave per- 
turbations of the most general type destabilise a static 
universe. We pointed out the link that can be forged 
between this homogeneous instability and the behaviour 
of the inhomogeneous gravitational wave spectrum by 
choosing modes with imaginary wave number. Our re- 
sults show that if the universe is in a neighbourhood of 
the Einstein static solution, it stays in that neighbour- 
hood, but the Einstein static is not an attractor (because 
the stability is neutral, with non-damped oscillations). 
Expansion away from the static state can be triggered 
by a fall in the pressure of the matter. Typically, expan- 
sion away from the static solution will lead to inflation. If 
inflation occurs, then perturbations about a Friedmann 
geometry will rapidly be driven to zero. The nonlinear 
effects (which will certainly be important in these models 
because of the initial infinite time scale envisaged) will 
be discussed in a further paper, as will other aspects of 
the spatially homogeneous anisotropic modes. 
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